Magnetoelectric Effect and Spontaneous Polarization in HoFe3(B03)4 and 

Hoo.5Ndo.5Fe3(B03)4 
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The thermodynamic, magnetic, dielectric, and magnetoelectric properties of HoFe3(B03)4 and 
Hoo.5Ndo.5Fe3(B03)4 are investigated. Both compounds show a second order Neel transition above 
30 K and a first order spin reorientation transition below 10 K. HoFe3(B03)4 develops a spontaneous 
electrical polarization below the Neel temperature (Tjv) which is diminished in external magnetic 
fields. No magnetic-field induced increase of the polarization could be observed in HoFe3 (603)4. In 
contrast, the solid solution Hoo.5Ndo.5Fe3(B03)4 exhibits both, a spontaneous polarization below 
Tiv and a positive magnetoelectric effect at higher fields that extends to high temperatures. The 
superposition of spontaneous polarization, induced by the internal magnetic field in the ordered 
state, and the magnetoelectric polarizations due to the external field results in a complex behavior 
of the total polarization measured as a function of temperature and field. 



PACS numbers: 75.30.Kz, 75.50.Ee, 75.80.+q, 77.70.+a 
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I. INTRODUCTION 

The rare earth iron borates, RFea (803)4 (R=rare 
earth, Y), belong to the trigonal system with space 
group R32.^ Their structure is similar to that of the 
mineral CaMg3(C03)4 (huntite), a trigonal trapezohe- 
dral structure that is one of the five trigonal types. 
This class of compounds has attracted recent attention 
because of a large magnetoelectric effect observed in 
GdFe3(B03)4^i^>ii^ and in NdFe3(B03)4.^ The noncen- 
trosymmetric structure makes these materials also an in- 
teresting candidate for optical applications based on their 
good luminescent and nonlinear optical properties i^i^ 
Thermodynamic and magnetic measurements have indi- 
cated a wealth of structural and magnetic phase transi- 
tions with ordering of Fe-spins as well as the rare earth 
moments involved. The majority of the RFcs (803)4 
compounds exhibit a structural phase transition from 
the high-temperature structure with R32 space group 
(No. 155) to the low-temperature P3i21 structure (No. 
152)fi2iiiii^ The structure change is accompanied by a 
distinct anomaly of the heat capacityJ^ii^ The R32 
P3i21 transition temperature, Tg, scales with the ionic 
radius of the rare earth and for some (R=Ho, Dy and 
Y) this transition takes place above room temperature. 
Other physical quantities, for example the dielectric con- 
stant, exhibit distinct and sharp anomalies at Tgi^J^ 

Magnetic measurements at high temperature reveal 
a strong antiferromagnetic (AFM) correlation between 
the Fe-spins resulting in an AFM transition at the Neel 
temperature, T^r, into the AFM2 phase below 40 K/^ 
Again, Tjv scales with the rare earth element decreas- 
ing from 40 K (Ho) to 22 K (La). Below Tat, the 
magnetic interaction of rare earth moments with the or- 
dered Fe-spins becomes stronger and results in a mag- 
netic polarization of the rare earth system and, in some 



rare cases, in a spin reorientation transition into the 
AFMl phase at lower temperatures that is driven by 
the magnetocrystalline anisotropy of the rare earth mo- 
ment. The spin reorientation transition was originally 
only observed in GdFe3 (803)4 at T5fl=9 K and its mag- 
netic field-temperature (H-T) phase diagram was exten- 
sively studiedi^ d'^d^d^ However, recent neutron scatter- 
ing experimentsii have found evidence for a similar spin 
reorientation in HoFe3 (803)4 at a slightly lower temper- 
ature of 5 K. The peculiarities of the crystallographic 
structure, the various possibilities for the Fe-spins to 
interact via direct or superexchange, and the coupling 
with the rare earth moments result in complex mag- 
netic structures and the cascade of phase transitions upon 
decreasing temperature, as observed in GdFe3(803)4ii^ 
The magnetic structure has been resolved for only a few 
rare earth iron borates through magnetic X-ray scatter- 
ing (R=Gd)i^ and neutron scattering (R=Nd, Tb, and 
Ho) experiments] "'^^d^'^" 

Since the discovery of a large magnetoelectric effect in 
gadolinium iron borate,^ the control of the electrical po- 
larization through magnetic fields in GdFe3 (803)4 has 
become the focus of interest. The complex magneto- 
electric H-T phase diagram was studied through polar- 
ization, dielectric constant, magnetic, and heat capacity 
mcasurementsi^'^'^ i^d'^d^i^i The combination of polariza- 
tion and magnetostriction experiments have proven the 
intimate correlation of the magnetoelastic and magneto- 
electric properties of the compound i^i^i^^ The sharp and 
sizable increase of the electrical polarization and the mag- 
netostriction effect is mainly observed at the spin reori- 
entation phase transition between the AFM2 and AFMl 
phases. The field-induced polarization (FIP) phase is 
well resolved in a plateau-like structure of the dielec- 
tric constant e{T) and the complete H-T phase diagram 
was resolved for different orientations of the external field 
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H.^ The magnetoelectric and magnetoelastic properties 
of GdFea (603)4 have been quahtatively described by a 
model taking into account the symmetry of the system 
and the lowest order expansion of the thermodynamic 
potential with respect to the magnetic exchange. 

The search for other rare earth iron borates with a sig- 
nificant magnetoelectric effect has been successful and 
NdFe3 (603)4 was discovered as the second candidate 
showing a large field-induced polarizationi^i^ Contrary 
to GdFca (603)4 no spin reorientation transition could be 
detected in NdFca (603)4 through magnetization, ther- 
modynamic, and neutron scattering experiments . ^^'^'^1^^ 
A small maximum of the magnetic susceptibility, ob- 
served at 6 K, was attributed to 3-dimensional AFM 
order of the Fe'^+ and Nd'^+ sublattices.^^ It is re- 
markable that the magnetic field-induced polarization 
in NdFe3 (603)4 exceeds the value in GdFe3 (603)4 
by nearly two orders of magnitude. The sizable 
magnetoelectric coupling found in GdFcs (603)4 and 
NdFe3 (603)4 are the motivation for the search for mag- 
netoelectric effects in other rare earth iron borates. Ho 
is an interesting candidate because of a relatively large 
magnetic moment of the rare earth ion. Strong interac- 
tions with the Fe-spins and possible magnetoelectric ef- 
fects are therefore expected in HoFe3 (603)4. In fact, the 
magnetic order in HoFe3 (603)4 below Tjv was recently 
explored and it was shown to involve long range order 
of both, Fe-spins and Ho-moments At lower tempera- 
tures, HoFcs (603)4 exhibits a spin reorientation similar 
to the one observed in GdFe3(603)4. However, investi- 
gations of the dielectric and magnetoelectric properties 
have not been reported yet. 

We have therefore synthesized large single crys- 
tals of HoFe3 (603)4 and of the solid solution 
Ndo.5Hoo.5Fe3 (603)4 and investigated their thermo- 
dynamic, magnetic, and magnetoelectric properties. 
We find that, even without external magnetic field, 
HoFe3(603)4 exhibits an electric polarization below Tjv 
which increases with decreasing temperature and sud- 
denly drops to zero at the spin reorientation transition. 
The magnetic field effect on the polarization and the H- 
T phase diagram is investigated. Ndo.5Hoo.5Fe3 (603)4 
shows a similar polarization effect upon decreasing 
temperature with an even higher value of the polar- 
ization as compared to HoFe3(603)4. In addition, 
Ndo.5Hoo.5Fe3 (603)4 also exhibits a large magnetoelec- 
tric effect in external fields along the a-axis with high 
values of the induced polarization. 

The paper is organized as follows: Section 2 provides 
the essential information about sample synthesis and ex- 
perimental equipment used for various measurements. 
The results are presented in Section 3 starting with 
data for HoFe3(603)4 (magnetic, thermodynamic, and 
dielectric data) and followed by corresponding data for 
Ndo.5Hoo.5Fe3 (603)4. Section 4 summarizes the main 
results and conclusions. 
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FIG. 1: (Color online) Magnetic susceptibility of 
HoFe3(B03)4 with the field oriented along the a-axis. 
The inset shows the derivative, d^a/dT, with the distinct 
anomaly at Tjv. 



II. EXPERIMENTAL 

Single crystals of HoFcs (603)4 and 
Ndo.5Hoo.5Fe3 (603)4 were grown using 6i2Mo30i2 
based fiux as described earlier 1^ The crystals were cut 
in different sizes and shapes to fit the demands of the 
various measurements. The orientation was determined 
by single crystal Laue X-ray diffractometry. For di- 
electric and polarization measurements thin plates (0.5 
to 1 mm thick) were prepared and electrical contacts 
were attached to two parallel faces using silver paint. 
For dielectric constant measurements the high-precision 
capacitance bridge AH2500A (Andeen-Hagerling) was 
employed. The electrical polarization was extracted 
by integrating the pyroelectric current measured by 
the K6517A electrometer (Keithley) upon variation of 
temperature or magnetic field. Temperature and field 
control were provided by the Physical Property Measure- 
ment System (PPMS, Quantum Design). The PPMS 
was also employed for heat capacity experiments. The 
magnetization was measured in a commercial SQUID 
magnetometer (MPMS, Quantum Design). Thermal 
expansion measurements were conducted using a home 
made high-precision capacitance dilatometer^i 



III. RESULTS AND DISCUSSION 

A. Magnetic phase diagram of HoFe3(B03)4 

Early investigations have shown that the structural 
transition from R32 to P3i21 in HoFe3 (603)4 happens 
well above room temperature. The critical temperatures, 
Ts ^ 420 K and T^ =i39 K, are among the highest 
of all rare earth iron borates.-iS The high magnetic or- 
dering temperature and the easy-plane character of the 
Fe'^"'' spin alignment was revealed also in optical measure- 
ments using Er'^"*' ions as a spectroscopic probe-^ The 
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FIG. 2: (Color online) Low-temperature magnetic susceptibil- 
ity of HoFe3(B03)4 with the field oriented along the c-axis. 



existence of a spin reorientation transition {Tsr —5 K) 
in HoFe3(B03)4, however, was detected only recently in 
neutron diffraction experiments at zero magnetic field.— 
The magnetic susceptibilities of HoFe3(B03)4, mea- 
sured perpendicular (xa) and parallel (xc) to the c-axis, 
clearly reveal the distinct anomalies at both magnetic 
transition temperatures (Figs. 1 and 2). The Neel tran- 
sition at T7v=38.5 K is reflected in a minute change of 
slope of the a-axis susceptibility which is clearly visible 
in the derivative, dxa/dT (inset in Fig. 1). Below 
the susceptibility shows a distinct dependence on the a- 
axis field. Ha. Xa increases with Hq below 5 kOe (Fig. 
I). This behavior is similar to the T-dependence of Xa of 
GdFe3(B03)4 which was interpreted as an effect of the 
in-plane ordering of the Fe-spins leading to a reduction of 
Xa between the Neel and spin reorientation temperatures 
at zero field. ^ At lower temperature, the spin reorienta- 
tion transition is marked by a sharp step of Xa- This step 
is positive with decreasing T for low magnetic fields Hq < 
1.2 kOe, however, it turns negative for Ha between 1.2 
kOe and 9 kOe. The sign change of this step-like anomaly 
is related to the magnetic order induced small decrease 
of Xa in the AFM2 phase since Xa is independent of the 
magnetic field in the AFMl phase (Fig. 1). The decrease 
of Xa at Tsr is different from the corresponding increase 
observed in GdFe3(B03)4. The origin of this difference 
may lie in the distinct magnetic structures of both com- 
pounds. The order of the Fe-spins is incommensurate 
with no long range order of the Gd- moments above T5 r 
in GdFe3(B03).fi^ but it is commensurate with the on- 
set of Ho-moment order at Tjv in HoFe3(B03)4.^^ Tsr 
quickly shifts to lower temperatures with increasing mag- 
netic field and the AFMl phase is completely suppressed 
for Ha > 10 kOe. The derived phase boundary between 
the AFM2 and AFMl phases is shown in Fig. 3 (bold 
symbols). 

The c-axis susceptibility shows no detectable anomaly 
at Tjv and very little field dependence in the AFM2 
phase. At the spin reorientation transition Xc suddenly 
drops by 60 % (Fig. 2). The resulting phase boundary is 



FIG. 3: (Color online) Low-temperature T-H phase diagram 
of HoFea (603)4 with the field aligned with the a-axis (bold 
symbols) and the c-axis (open symbols). The lines serve as a 
guide to the eye. 



displayed in Fig. 3 (open symbols). The large decrease 
of Xc in the AFMl phase reflects the increased correla- 
tion and stiffness of the Fe-spins after their rotation into 
a collinear spin structure oriented along the c-axis that is 
induced by the strong coupling to the Ho moments with 
a large uniaxial magnetic anisotropy. The magnetic sus- 
ceptibility data reveal a minute temperature hysteresis 
of about 0.1 K at Tsr indicating the first order nature 
of the phase transition from the AFM2 to the AFMl 
phase. The first order nature of the spin reorientation 
phase transition was also suggested from the step-like 
change of the Ho and Fe moments derived from powder 
neutron scattering experiments The phase diagram is 
further explored by isothermal M-H measurements show- 
ing a metamagnetic transition (sudden increase of the 
magnetization M(H) with field) at the spin reorienta- 
tion transition. The critical fields are consistent with 
the phase boundaries of Fig. 3 derived from temperature 
dependent measurements. Similar to the x(T) measure- 
ments the M(H) data show a small field hysteresis across 
the AFM2^AFM1 phase boundary confirming the first 
order character of the transition. 

The thermodynamic signature of the two magnetic 
phase transitions is revealed in heat capacity measure- 
ments. Figure 4 shows the heat capacity of HoFe3(B03)4 
together with Cp of the related compounds GdFe3(B03)4 
and ErFe3(B03)4 in a wide temperature range. The in- 
set displays the low-temperature Cp of HoFe3(B03)4 near 
the spin reorientation transition on an expanded temper- 
ature scale. The Neel transition is characterized by a A- 
shaped pronounced peak of Cp(T) at 38.5 K suggesting 
the transition into the AFM2 phase to be a second or- 
der phase transition. This is consistent with the steady 
increase of the Fe sublattice magnetization below T^r ob- 
served in neutron scattering experiments. However, the 
derived Fe magnetic moment does not follow closely the 
expected Brillouin function which was attributed to the 
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FIG. 4: (Color online) Comparison of heat capacities of 
HoFe3(B03)4, GdFe3(B03)4, and ErFe3(B03)4. The inset 
shows the low-temperature heat capacity of HoFe3(B03)4 
near Tsr- 



coupling of the Fe-spins with the Ho-moments leading to 
a simultaneous order of both magnetic subsystems at Tat. 
The low-temperature peak of Cp at Tsr is extremely 
narrow as is expected at a first order transition. Be- 
low Tat, but above Tsr, Cp/T of HoFe3(B03)4 shows 
an enhancement with respect to similar data for other 
rare earth iron borates, e.g. ErFcs (803)4 also shown 
in Fig. 4. This enhancement reflects the strong coupling 
between the rare earth moments and the Fe spin order re- 
sulting in the spin reorientation transition and the order 
of the rare earth moments. An equivalent enhancement 
of Cp/T was also observed in GdFea (603)4 (Fig. 4), the 
second RFe3 (603)4 with a spin reorientation transition.^ 
The increased eirtropy associated with the enhanced heat 
capacity at low T was attributed to the eight-fold degen- 
eracy of the ground state level of the Gd'^^ ion in the 
corresponding crystal fieldf^^ The heat capacity data for 
all three compounds are compared in Fig. 4. At the 
lowest temperatures a Schottky-type contribution to the 
heat capacity is most pronounced in ErFe3(B03)4 and 
GdFe3 (603)4, but barely visible in HoFe3(B03)4. 



B. Magnetoelectric effect and polarization in 

HoFe3(B03)4 

1. Dielectric constant m magnetic fields 

Dielectric anomalies below the Neel transition have 
been observed in GdFcs (603)41^ A small but distinct in- 
crease of the a- axis dielectric constant, Eq, and its sudden 
drop at the spin reorientation transition was explained by 
the coupling of the magnetic order to the lattice. While 
the peak-like increase of £a in GdFe3 (603)4 is less than 
1 % of its value at Tjv, the a-axis dielectric constant 
of HoFe3 (603)4 (Fig. 5) shows a colossal enhancement 
of nearly 100 % below Tm- This huge increase from 



ea=20 to £0=37 indicates that the spin-lattice coupling in 
HoFe3(603)4 is two orders of magnitude stronger than 
the spin-lattice interaction in GdFe3(603)4. With de- 
creasing temperature Sa passes through a maximum and 
drops suddenly back to its high-tcmperature value of 20 
at the spin reorientation transition temperature, Tsr, as 
shown in Fig. 5. 

The peak of Sa is quickly suppressed by magnetic fields 
oriented along the a-axis (Fig. 5a). At 4 kOe, the peak 
height of Ea has dropped by 90 % and at 10 kOe it is com- 
pletely absent. This extraordinary magneto-dielectric ef- 
fect shows the intricate correlation between the magnetic 
order and the dielectric properties of HoFe3 (603)4. With 
the magnetic field oriented along the c-axis (Fig. 5b) the 
rise of £a below Tat persists to higher field values, up to 
He ~8 kOe, beyond the critical field above which the spin 
reorientation transition is suppressed (see Fig. 3). The 
maximum Ea —40 at 8 kOe is twice as large as its value 
near or above T^v- With further increasing field, how- 
ever, the maximum of Ea decreases quickly but a small 
enhancement below Tat is still detectable at c-axis fields 
as high as 40 kOe, as shown in Fig. 5b. It is remarkable 
that the main increase of Ea in the AFM2 phase hap- 
pens well below the Neel temperature at about 20 K and 
that Ea suddenly decreases at Tsr assuming a value in 
the AFMl phase that is close to Ea of the paramagnetic 
phase above Tat. This shows that the coupling between 
the magnetic order and the lattice is most significant in 
the AFM2 phase at zero magnetic field but it is negli- 
gibly small in the AFMl phase. The c-axis dielectric 
constant, Ec, was also measured and it shows a similar 
enhancement below Tat, however, the magnitude of the 
peak is much smaller and it amounts to only 3 % of its 
base value. 

A similar large magneto-dielectric effect has been ob- 
served in some rare earth manganites, for example in 
orthorhombic HoMnOs below the Neel temperature.— It 
was later shown that the sharp increase of e was asso- 
ciated with the development of a spontaneous electric 
polarization.— Therefore, we decided to search for pyro- 
electric and magnetoelectric effects below the Neel tem- 
perature of HoFe3(603)4. 



2. Polarization at zero magnetic field 

The pyroelectric current measured between two paral- 
lel electrodes attached to a sample is proportional to the 
change of the sample's polarization due to a change of 
temperature or magnetic field. In general, a bias voltage 
can be applied during the measurement or it can be used 
to align domains in a ferroelectric state. The integration 
of the pyroelectric current provides a quantitative mea- 
sure of the polarization change. However, special pre- 
cautions have to be taken to ensure that the measured 
current actually reflects the polarization change and is 
free from artifacts such as electrical transport currents 
or capacitive contributions due to the sample's dielectric 
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FIG. 5: (Color online) Dielectric constant, Sa, of 
HoFea (603)4 at different magnetic fields H||a (a) and H||c 
(b). 



properties. The current measured in an experiment can 
be expressed by equation (1). 



1 = — + C— + —V + A— 
R dt dt dt 



(1) 



V is the apphed bias voltage, R and C are the resis- 
tance and capacitance of the sample, respectively, and P 
is the intrinsic electrical polarization. A is the contact 
area. The first term is the resistive current that can be 
avoided if the sample has a very high resistance or the 
measurement is conducted at zero bias, V=0. The effects 
of charge trapping and release with changing temperature 
usually are not significant for insulating materials. The 
second term is related to the charging current of a ca- 
pacitor and will be zero if the applied voltage is constant 
with time. The third contribution in (1) may cause a cur- 
rent if the dielectric constant (capacitance) of the sample 
changes significantly with temperature (as for example Sa 
in Fig. 5) and a bias voltage is applied. The resulting 
current can be estimated from the dielectric data and is 
negligible in most experiments. It can be avoided com- 
pletely if no voltage bias is applied. The last term of (1) is 
the pyroelectric current due to the change of the sample's 
polarization. In the current experiments all polarization 
measurements have been conducted at zero bias voltage, 



V=0, which ensures that the first three terms of (1) do 
not cause any additional contribution to the current so 
that the integration of the measured pyroelectric current 
yields the sample's polarization change. Measurements 
are commonly conducted by cooling the sample with an 
applied electrical bias voltage and measuring the pyro- 
electric current in zero bias upon warming. 

In HoFca (603)4 we found that even in zero- voltage 
cooling the crystal developed a pyroelectric current aris- 
ing below Tjv at about 25 K. In passing through the spin 
reorientation transition the current shows a sharp nega- 
tive peak indicating a sudden drop of the polarization. 
Integrating the current reveals the associated change of 
polarization shown in Fig. 6 ((a): a-axis and (b): c-axis 
polarization) for cooling and warming. Upon decreasing 
temperature, the polarizations measured along the a- and 
c-axes gradually increase in the AFM2 phase and reach 
a maximum of 60 /iC/m^ and 90 /iC/m^, respectively, at 
TsR- A sharp drop of Pa,c to a small value (presumably 
zero, the finite value shown in Fig. 6 is within the error 
limits of the measurement and integration) indicates the 
transition into the AFMl phase. The maximum values of 
Pa and Pc are comparable with the magnitude of the po- 
larization in the ferroelectric phase of other multiferroics, 
for example Ni3V208 and MnW04i^i2^ However, there 
is a major difference in the temperature dependence of 
the polarization of HoFe3(B03)4 in comparison with the 
typical behavior observed in the majority of multiferroic 
ferroelectric compounds. In the latter compounds the 
ferroelectric polarization quickly increases below a mag- 
netic phase transition and, in most cases, the direction 
of P is determined by the direction of the bias voltage 
applied during electric-field cooling. In the present data 
for IIoFe3(B03)4 there is no significant change of P(T) 
right below the magnetic transition temperature (38.5 
K) but Pq and Pc smoothly increase below a much lower 
temperature (~ 25 K). Applying a pohng voltage (± 150 
V) had no sizable effect on P(T) and the polarization 
could not be reversed. Pyroelectric current data mea- 
sured upon cooling with a poling voltage applied as well 
as upon warming after the bias voltage was reduced to 
zero at 6 K are almost identical to the zero-bias data 
shown in Fig. 6. This indicates a strong preference of 
the direction of Pa and Pc which may be determined by 
structural or magnetic domains in contrast to most mul- 
tiferroic ferroelectric systems. A clear thermal hysteresis 
of P(T) between 15 K and 25 K (indicated by arrows 
in Fig. 6) provides further evidence that the value of 
the polarization in the AFM2 phase is sensitive to the 
magnetic domain structure. 

The most dramatic change of the dielectric constant (~ 
100 % increase. Fig. 5) was observed along the a-axis. 
The maximum of the polarization Pa(T) is also larger by 
about 50 % than the maximum of Pc (Fig. 6) and the 
thermal hysteresis is even more significant and extends 
between Tg/j and about 30 K. To investigate the effects 
of electric fields (bias voltage) on the polarization Pa the 
pyroelectric current was measured upon cooling to 6 K 
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FIG. 6: (Color online) Polarization of HoFea (603)4 measured 
along (a): the a-axis and (b): the c-axis in zero electric and 
magnetic fields. The arrows distinguish the cooling and warm- 
ing data. 



(above the spin reorientation transition) with a bias of ± 
200 V apphed. At 6 K the external voltage was set to zero 
and, after short-circuiting the contacts, the remaining 
polarization was measured upon warming from 6 K to 
50 K. The data shown in Fig. 7 reveal a large apparent 
change of the polarization due to the external voltage. 
At + 200 V cooling. Pa (+200 V) at 6 K is enhanced by 
50 % with respect to the zero-bias Pa(0 V) and at - 200 V 
cooling, Pa(-200 V) is reduced to a fraction of less than 
20 %. The solid lines in Fig. 7 represent the polarization 
data obtained from integrating the pyroelectric current 
upon cooling in + 200 V (top hne), V (center line), 
and - 200 V (bottom line). After the release of the bias 
voltage at 6 K, however, the polarization assumes the 
zero bias value close to 70 /zC/m^, as indicated by the 
vertical dotted arrows. Upon warming Pa follows closely 
the zero-bias temperature dependence (dashed lines in 
Fig. 7). The large change of Pa at 6 K upon release of 
the bias voltage seems to indicate a significant tunability 
of the polarization in electric fields. 

However, the dielectric constant Sa also increases sig- 
nificantly below T^v (Fig. 5) and the current contribution 
from the third term of equation (1) has to be taken into 
account in measurements conducted in an external bias 
electric field to determine the intrinsic polarization val- 
ues. From the dielectric measurements we determine the 
capacitance increase in the AFM2 phase as about 1 pF 
which results in an additional charge of ± 2*10~^" C due 
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FIG. 7; (Color online) a-axis polarization of HoFe3(B03)4 
measured upon cooling (solid lines) with applied bias voltage 
of -1-200 V, V, and - 200 V. Data upon warming (dashed 
lines) are collected after release of the bias voltage at 6 K. 



to a bias voltage of ± 200 V. The change of the appar- 
ent polarization at 6 K in Fig. 7 (indicated by vertical 
dotted arrows) amounts to 2.2*10-1'' C and -2.7*10-i" 
C for a bias voltage of -I- 200 V and - 200 V, respectively. 
After correcting for these additional contributions to the 
pyroelectric current, the values for Pa appear to change 
only very little in an electric field. Therefore, we con- 
clude that for the a-axis pyroelectric measurement the 
zero-bias data (dashed lines in Fig. 7) are representa- 
tive for the intrinsic polarization, Pa(T). All subsequent 
pyroelectric measurements in magnetic fields have been 
conducted in zero electric field to avoid the additional 
contribution from the third term in (1) to the pyroelec- 
tric current measurement. It should be noted that this 
problem basically does not exist in the c-axis measure- 
ments (Fig. 6a) since the increase of Sc in the AFM2 
phase is negligibly small, only 3 % of its value near T^r. 

The appearance of a spontaneous polarization P(T) in 
the AFM2 phase and its relative insensitivity to elec- 
tric fields in HoFea (603)4 is unique among the rare 
earth iron borates. Changes of electrical polarization 
in RFe3 (603)4 have previously only been observed with 
the application and change of an external magnetic field 
for R=Gd^' and R=NdF^ The spontaneous polarization 
in HoFe3(B03)4 has its possible origin in the coupling to 
the internal magnetic field that is associated with the or- 
der of the iron spins as well as the Ho moments. It should 
be noted that a polarized state is not forbidden by the 
crystal symmetry because the P3i21 space group is non- 
centrosymmetric. The magnetic order below Tjv and its 
strong coupling to the lattice results in the development 
of a macroscopic electrical polarization. The preferred 
orientation of this polarization and the observation that 
it cannot be reversed by electric fields is a result of the 
non-centrosymmetric structure that allows for one direc- 
tion of the polarization only within a structural and/or 
magnetic domain. Magnetic domains in the AFM2 phase, 
however, also play an important role as evidenced by the 
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thermal hysteresis observed in the polarization data (Fig. 
6) and it lends further support to the role of the inter- 
nal magnetic field/order and its mutual interaction with 
the polarization. From the measured polarization values 
along the a- and c-axes we conclude that the intrinsic 
polarization direction is between both major crystallo- 
graphic orientations, at an angle of about 55 ° with the 
c-axis. The origin of the intermediate orientation must be 
sought in the peculiarities of the magnetic structure. Ac- 
cording to recent neutron scattering experiments^^- there 
exist two inequivalent Fe-sites. In the polarized AFM2 
phase the Fe(l)-spins and the Ho-moments are aligned 
with the a-axis, however, the spins on the Fe(2) sublat- 
tice are tilted towards the c-axis. The estimated val- 
ues of the X- and z-components of the Fe(2)-spins at 5 
K are 4 /i^ and 2.5 fiB, respectively. The angle of the 
Fe(2)-spins with the c-axis is therefore approximately 58 
°, close to the polarization orientation discussed above. 
While it seems unlikely that this coincidence is acciden- 
tal it rather demonstrates the intricate coupling of the 
intrinsic polarization with the specifics of the magnetic 
order, supposedly of the Fe(2) sublattice. 

The sharp drop of Pq and Pc at the spin reorientation 
transition indicates that the observed polarization is only 
allowed in the AFM2 phase with most of the magnetic 
moments oriented along the a-axis (some iron spins have 
a c-axis component, however).i- In the low-temperature 
AFMl phase all iron spins and half of the Ho moments 
are aligned with the c-axis and the remaining Ho mo- 
ments exhibit a tilt toward the a-axis. The particular 
magnetic order apparently does not couple to the lattice 
to induce a macroscopic polarization. The microscopic 
details and the symmetry constraints have yet to be in- 
vestigated. Forthcoming studies may also include the in- 
vestigation of the dielectric relaxation behavior to gather 
complimentary knowledge about the nature of the polar- 
ized statci^ 



3. Magnetic field effects on the polarization 

To investigate the effects of external magnetic fields 
on the polarization, pyroelectric measurements have been 
conducted in zero electric field upon warming from the 
lowest temperatures (2 K) in a constant magnetic field H 
applied along the a- and c-axes. Results for the electrical 
polarization, Pc, are shown in Fig. 8. With H oriented 
along the a-axis Pc is quickly diminished in even moder- 
ate magnetic fields (Fig. 8a). The sharp drop of Pc to 
zero at Tsr persists up to the critical field above which 
the AFMl phase is completely suppressed. At higher 
fields (10 kOe in Fig. 8a) Pc remains finite to the low- 
est temperatures although the overall magnitude is only 
about 15 % of its maximum value at zero field. With 
the magnetic field applied along the c-axis (Fig. 8b) the 
spin reorientation transition is quickly suppressed at low 
fields but the polarization retains its high values up to 
8 kOe. At this field Pc remains large to the lowest tem- 




T(K) 

FIG. 8: (Color online) c-axis polarization of HoFea (603)4 at 
different magnetic fields H. (a): H |[ a-axis and (b): H || c-axis. 



peratures. Only further increasing field reduces the mag- 
nitude of the polarization but at a lower rate as for the 
a-axis field. At 30 kOe the low-T polarization is still as 
large as 8 /iC/m^. The phase boundary separating the 
AFM2 and AFMl phases as derived from the polariza- 
tion drop (Fig. 8) is in perfect agreement with the phase 
diagrams shown in Fig. 3. 

The polarization measured along the a-axis, Pci(T), is 
shown for different magnetic fields oriented along a and 
c in Fig. 9. Pa(T) reveals a similar suppression in exter- 
nal magnetic fields as Pc(T) discussed above. The a-axis 
magnetic field has the strongest effect on reducing P^ 
and even 2 kOe are sufficient to decrease Pa by 50 %. 
The c-axis magnetic field suppresses the spin reorienta- 
tion transition first while maintaining a relatively high 
magnitude of Pa up to about 10 kOe. Only higher fields 
He result in a slow decrease of Pa which reaches about 
10 % of the zero-field value at Hc=40 kOe. The hump 
of Pa that develops in magnetic fields between 15 K and 
20 K (Fig. 9) is associated with the temperature hys- 
teresis in this range as shown in the zero-H polarization 
in Fig. 6. In magnetic fields the hysteretic region (ex- 
tending from Tsr to about 25 K at H=0, Fig. 6a) shifts 
toward higher temperature. For better clarity the po- 
larization measured upon cooling is not included in Fig. 
9. It seems surprising that the intrinsic zero-field polar- 
ization of HoFca (603)4 is quickly suppressed by external 
fields whereas in the related compound NdFca (603)4 the 
magnetic field effect is significant and results in a large 
increase of the polarization. The major differences lie 
in the magnetic order of rare earth and iron moments 
below the Neel temperature. In HoFca (603)4 the spins 
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T(K) 

FIG. 9: (Color online) a-axis polarization of HoFe3(B03)4 at 
different magnetic fields H. (a): H || a-axis and (b): H || c-axis. 

on the Fe(l) sublattice are coUinear (oriented along the 
a-axis) with the Ho-momcnts but noncoUinear with the 
remaining spins on the Fe(2) sublatticejii In contrast, 
all Fe-spins in NdFea (603)4 are collinear with an easy 
plane anisotropy and the Nd-moments are noncoUinear 
with the Fe-spins^i2, These differences in the magnetic 
structures may be responsible for the different magneto- 
electric properties. 

4- Thermal expansion anomalies 

The magnetoelectric coupling in multiferroic and mag- 
netoelectric compounds are associated with ionic dis- 
placements at phase transitions that commonly result in 
macroscopic distortions or lattice strain along the crys- 
tallographic axeSi ^^'^^i'^^'^'^i'^^ The relative change of the 
lattice parameters and the volume of HoFea (603)4 was 
therefore measured and is shown in Fig. 10. The refer- 
ence temperature for both axes was chosen at 55 K. The 
dashed vertical lines indicate the two magnetic transi- 
tions. At T^r the temperature dependence of the a- and 
c-axes both exhibit a minute change of slope, as is ex- 
pected for a second order magnetic phase transition. At 
Tgfl, however, a sharp increase of c and the drop of a 
are characteristic for a first order phase transition. The 
significant changes of the lattice parameters at Tsr are 
consistent with the abrupt disappearance of the electri- 
cal polarization and the sharp changes of the magnetic 
susceptibilities (Figs. 1 and 2). The large anomalies of 
the lattice parameters reveal the strong magnetoelastic 
effect in HoFea (603)4. It is interesting that the thermal 
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FIG. 10: (Color online) Thermal expansion of the a- and c- 
axes and the volume of HoFe3(B03)4. The two magnetic 
phase transitions are marked by vertical dashed lines. 

expansivities are strongly anisotropic and the c-axis ex- 
pansivity is negative over most of the temperature range 
shown in Fig. 10. Similar negative expansivities have 
been observed, for example, in multiferroic HoMnOg^ 
and DyMn205.34 It was understood as a signature of 
strong magnetic correlations and spin-lattice interactions 
in connection with large magnetic anisotropy. The c-axis 
length exhibits a maximum at about 15 K. At lower tem- 
peratures both, c and a, decrease faster with decreasing 
temperature approaching the instability of the magnetic 
structure at the spin reorientation phase transition. The 
relative volume change (Fig. 10) clearly shows a rapid 
decrease below 15 K leading to the magnetic instability 
at TsR=i.6 K. 

C. Magnetic phase diagram, magnetoelectric 
effect, and polarization in Hoo.5Ndo.5Fe3(B03)4 

NdFcs (603)4 does not show any evidence for a spin 
reorientation phase transition at low temperatures^i^ 
but it exhibits a large magnetoelectric effect<^ The solid 
solution of NdFe3 (603)4 and IIoFe3 (603)4 is therefore 
of interest with regard to the magnetic phase diagram 
and the magnetoelectric properties. 



1. Magnetic phase diagram 

The magnetic susceptibility was measured in 
Hoo.5Ndo.5Fe3 (603)4 for both field orientations, 
Hq and He. Xa and Xc show a sizable step-like change 
E^t Tsfl ~9 K, clearly revealing the existence of a spin 
reorientation transition. While the drop of Xc at Tsr 
is large in passing into the AFMl phase and it amounts 
up to 80 %, the change of Xa is moderate and less 
than 20 %, similar to HoFe3(603)4. However, the 
anomaly of Xa changes sign at about 2 kOe. 6elow 2 
kOe, upon decreasing temperature, the step of Xa is 
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positive whereas above 2 kOe it is negative. At 2 kOe 
the anomaly of Xa at Tsr is a change of slope with no 
sudden change, as shown in Fig. 11a. A similar behavior 
is also observed in the polarization in magnetic fields 
applied along the a-axis (next section). It is remarkable 
that Tsr of Hoo.5Ndo.5Fe3 (603)4 is twice as large 
as that of HoFes (603)4. However, it is comparable 
with Tsr of GdFe3(B03)4.- The larger value of Tsr 
indicates that Nd stabilizes the low temperature AFMl 
phase and Tsr passes through a maximum in the phase 
diagram of the soHd solution Hoi-ajNda; Fes (603)4- 
The AFMl phase in Hoo.5Ndo.5Fe3 (603)4 is also more 
stable with respect to magnetic fields as compared to 
H0FC3 (603)4. The magnetic phase diagram derived 
from a- and c-axis magnetization measurements is shown 
in Fig. 12. The phase diagram includes data of x(T) 
acquired upon cooling and warming and results from 
isothermal magnetization measurements with increasing 
and decreasing fields. In comparison to HoFe3(603)4 
the AFMl phase of Hoo.5Ndo.5Fe3 (603)4 is more stable 
with respect to temperature as well as magnetic field. 
The zero-temperature critical fields of the AFMl — + 
AFM2 transition are estimated as 22 kOe and 13 kOe 
for H II a and H || c, respectively. 

The thermodynamic signatures of both magnetic phase 
transitions, paramagnetic — > AFM2 and AFM2 
AFMl, are resolved in distinct anomalies of the heat ca- 



FIG. 13: Heat capacity of Hoo.5Ndo.5Fe3(B03)4. 

pacity shown in Fig. 13. The A-shaped peak at Tjv— 32 K 
is characteristic for a second order phase transition. The 
first order spin reorientation transition is clearly marked 
by the sharp peak at about 9 K. While the peak at Tsr 
appears on top of a broad shoulder there is no clear max- 
imum of Cj,/T at lower temperatures that could indi- 
cate a Schottky type anomaly as observed for example 
in NdFe3(603)4^^ or GdFe3(603)4 (Fig. 4). The major 
difference to HoFe3(603)4 discussed above is the signifi- 
cantly lower Neel temperature of the iron spin order and 
the larger spin reorientation temperature. 

2. Magnetoelectric properties of Hoo.5Ndo.5Fe3 (303)4 

At zero magnetic field the spontaneous polarization 
measured along the a-axis, Pa(T), arises below 40 K in 
Hoo.5Ndo.5Fe3(603)4 and it reaches a maximum at the 
spin reorientation transition (Fig. 14). The maximum 
value of 120 ^C/rn^ exceeds the corresponding value of 
HoFe3(603)4 by a factor of two. The transition into 
the AFMl phase is accompanied by a sharp drop of Pa, 
however, P^ remains finite (about 40 /iC/m^) in the low 
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temperature AFMl phase. This is in distinct contrast to 
the behavior in HoFca (603)4 (Fig. 9), where Pa drops 
to zero in the AFMl phase. With the external magnetic 
field applied along the c-axis (Fig. 14b), the spin reorien- 
tation transition temperature is reduced resulting in an 
enhancement of the maximum of Pa(T) to 160 iiC/m? 
for fields up to 20 kOe. At higher magnetic fields He 
the polarization decreases continuously over the whole 
temperature range and its maximum value drops to 20 
IxC/m^ at Hc=7 T. 

The magnetic field dependence of Pa of 
Hoo.5Ndo.5Fe3 (603)4 for H||c is qualitatively simi- 
lar to the field effect in HoFe3 (603)4 as is also obvious 
from the isothermal polarization measurements shown 
in Fig. 15b. The polarization Pa(Hc) of Fig. 15 was 
derived from current measurements at constant temper- 
ature and increasing field. The value of Pq at zero field 
was determined from the T-dependent data (Fig. 14) at 
Hc=0 as the reference value. At 5 K the polarization 
Pa is finite and nearly constant for small fields but it 
increases sharply at 10 kOe because of the transition 
from the AFMl to the AFM2 phase. With further 
increasing field Pa (5 K) decreases quickly, consistent 
with the T-dependent data of Fig. 14b. At higher 
temperatures (> 10 K) Pa starts higher at zero field 
because of its larger value in the AFM2 phase. However, 
the external magnetic field H^ reduces Pa quickly as 
shown in Fig. 15b. The critical temperatures and fields 
of the spin reorientation phase transition determined by 
the sharp change of the polarization (Figs. 14b and 15b) 
coincide with the c-axis magnetic phase diagram of Fig. 
12 derived from magnetization measurements. 

The effect of the c-axis magnetic field on the polar- 
ization of Hoo.5Ndo.5Fe3 (603)4 is similar to the field- 
induced suppression of the spontaneous polarization of 
H0FC3 (603)4. However, the picture changes completely 
if the magnetic field is applied along the hexagonal a- 
axis. The polarization Pa(T) at various Ha is shown 
in Fig. 14a. In zero magnetic field Pa(T) arises below 
Tat, exhibits the maximum at Tsr and drops sharply 
to a smaller value in the AFMl phase, as discussed 
above. With increasing magnetic field, however, the 
positive Pa(T) is suppressed, it becomes negative, and 
passes through a minimum with decreasing temperature, 
as shown in Fig. 14a for Ha=10 kOe. The sudden de- 
crease of Pa(H=0) at T^sR is inverted in fields above 2 
kOe and Pa increases sharply with the transition into the 
AFMl phase. This behavior of Pa at ^SR is consistent 
with the anomaly of the magnetization discussed in the 
previous section demonstrating the close relation of the 
spontaneous polarization and the magnetization. 

With the magnetic field increasing further above 10 
kOe the polarization Pa (T) starts increasing again in the 
whole temperature range even above the AFM Neel tem- 
perature. Pa is all positive for Ha >20 kOe and it in- 
creases to large values, Pa >800 /xC/m^, at 70 kOe and 
2 K (Fig. 14a). This high-field behavior is similar to the 
magnetoelectric effect observed in NdFe3 (603)4 and can 
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FIG. 14; Temperature dependence of Pa(T) of 
Hoo.5Ndo.5Fe3(B03)4 in different magnetic fields, (a); 
H II a-axis and (b): H || c-axis. 
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FIG. 15: Magnetic field dependence of the polarization Pa at 
difi'erent temperatures, (a): H || a-axis and (b): H || c-ax;is 
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be attributed to a quadratic magnetoelectric coupling as 
discussed recently 

The complex temperature and field dependence of 
Pa(T) shown in Fig. 14a is also reflected in the isother- 
mal polarization data, Pa(Ha). The results at different 
temperatures are displayed in Fig. 15a. At zero and 
small field Pa is positive and determined by the sponta- 
neous polarization as discussed above. With increasing 
field the spontaneous Pa is suppressed, similar to the ef- 
fect of a c-axis field (Figs. 14a and 15a) and the field ef- 
fect in HoFe3 (603)4 (Section B). The magnetoelectric ef- 
fect, since quadratic in Ha,'' is negligibly small. The step- 
like anomaly of Pa near 22 kOe (2 K and 5 K data of Fig. 
15a) is consistent with the temperature anomaly of Pa 
at TsR (Fig. 14a). However, at higher field the magne- 
toelectric coupling dominates, similar to NdFe3(B03)4. 
The magnetoelectric effect at first decreases the polar- 
ization and Pa (Ha) changes sign and becomes negative 
above about 4 kOe. With further increasing Ha the mag- 
netoelectric effect changes sign and results in another sign 
reversal of Pa(Ha) at about 20 kOe. This effect was also 
observed in NdFe3 (603)4 and it was explained by the 
competition of the external magnetic field with the ex- 
change field;^ The major difference of the polarization 
behavior of NdFe3 (603)4 and Hoo.5Ndo.5Fe3 (603)4 is 
the existence of a polarization Pa in zero magnetic field 
which is aligned with the high-field magnetoelectric po- 
larization. The spontaneous polarization and the field 
induced magnetoelectric polarization are superimposed 
and result in a twofold sign reversal of the total polariza- 
tion upon increasing external magnetic field between 10 
K and 20 K, as shown in Fig. 15a. 

We have also conducted pyroelectric measurements 
along the hexagonal c-axis. The detected pyroelectric 
current was very small and any possible component of 
the polarization in zero and high magnetic fields was be- 
yond the resolution limits of the experimental procedure. 



IV. SUMMARY AND CONCLUSIONS 

We have investigated the thermodynamic, mag- 
netic, dielectric, and magnetoelectric properties 
of HoFe3(603)4 and Hoo.5Ndo.5Fe3(603)4. Both 
compounds show an antiferromagnetic Neel transi- 
tion at TAr=38.5 K (HoFe3(603)4) and Tjv=32 K 
(Hoo.5Ndo.5Fe3 (603)4) into the AFM2 phase and 
a spin reorientation phase transition at Tsr into a 
low-temperature AFMl phase. The AFMl phase is 
suppressed in external magnetic fields oriented along the 
a- and c-axes. The thermodynamic signature of both 
magnetic transitions is shown in distinct anomalies of the 
heat capacity indicating a second order transition at Tn 
and a first order phase transition at Tsr- The magnetic 
phase diagram of both compounds is completely resolved 
up to 70 kOe. 

The dielectric and magnetoelectric experiments on 
HoFe3 (603)4 reveal giant magneto-dielectric effect and 



the existence of a spontaneous polarization below the 
Neel temperature which suddenly drops to zero at Tsr- 
This effect has not been observed before in rare earth 
iron borates. The polarization in HoFes (603)4 has com- 
ponents along the a- and c-axes indicating its inter- 
nal orientation in between both crystallographic orienta- 
tions. This polarization only arises in the AFM2 phase of 
HoFcs (603)4 which shows that P is tied to the magnetic 
order and the resulting internal magnetic field of this 
phase. In contrast to NdFe3 (603)4 external magnetic 
fields do mainly suppress the polarization and the magne- 
toelectric effect is negative, dP/dH<0. A symmetry anal- 
ysis similar to NdFe3 (603)4 (Refj^) may help to clarify 
the origin of the zero-field polarization of HoFe3 (603)4. 
This analysis has to take into account the specific mag- 
netic order of both, Fe and Ho subsystems, as known 
from neutron scattering data.— 



Substituting 50 % of Ho with Nd results in a reduc- 
tion of Tjv to 32 K and an increase of the stability 
range of the AFMl phase. The spontaneous polariza- 
tion at zero magnetic field is also observed, mainly along 
the a-axis, with an increase of the maximum polariza- 
tion (at zero magnetic field) by 50 % as compared to 
the values in HoFcs (603)4. Similar to HoFe3 (603)4, 
magnetic fields along the c-axis suppress the polariza- 
tion of Hoo.5Ndo.5Fe3 (603)4. However, the a-axis mag- 
netic field couples to the polarization and the magneto- 
electric effect induces large values of the polarization Pa 
(up to 900 /iC/m^ at 70 kOe). The temperature and 
field dependence of Pa is very complex with several sign 
reversals of the polarization and step-like anomalies at 
the spin reorientation phase transition. The detailed fea- 
tures can be understood as a superposition of the spon- 
taneous polarization induced by internal fields and the 
magnetoelectric effect due to the external magnetic field. 
The solid solution Hoo.5Ndo.5Fe3 (603)4 combines fea- 
tures observed in both parent compounds, HoFcs (603)4 
and NdFe3 (603)4. 



Further investigations of the magnetic structure, the 
order of Fe-spins, Ho- and Nd-moments, and the mag- 
netoelectric interactions are necessary to gather a ba- 
sic understanding of the complex physical properties 
of Hoo.5Ndo.5Fe3 (603)4 and related compounds. The 
present investigation may be extended to include the 
solid solution of other members of the rare earth iron 
borate compounds. One interesting candidate could be 
the solid solution of Hoi_j;Gd2;Fe3(603)4. GdFe3(603)4 
is the only second compound in the rare earth iron bo- 
rate system to show a significant magnetoelectric ef- 
fect and, therefore, the magnetoelectric properties of 
Hoi_a;Gda;Fe3 (603)4 should be interesting too. A thor- 
ough investigation of Hoi_a:Gd2;Fe3(603)4 could even 
help to arrive at a deeper and more comprehensive under- 
standing of the current results on Hoo.5Ndo.5Fe3(603)4. 
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